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We study the quantum phase transition of an itinerant antiferromagnet with cubic anisotropy in the 
presence of quenched disorder, paying particular attention to the locally ordered spatial regions that 
form in the Griffiths region. We derive an effective action where these rare regions are described in 
terms of static annealed disorder. A one loop renormalization group analysis of the effective action 
shows that for order parameter dimensions p < 4 the rare regions destroy the conventional critical 
behavior. For order parameter dimensions p > 4 the critical behavior is not influenced by the rare 
regions, it is described by the conventional dirty cubic fixed point. We also discuss the influence of 
the rare regions on the fluctuation-driven first-order transition in this system. 



I. INTRODUCTION 

Quenched disorder can have very drastic influences on 
the critical behavior of a system undergoing a continuous 
phase transition. According to the Harris criterion^ the 
critical behavior of a clean system is unaltered by dis- 
order, if the correlation length critical exponent v obeys 
the inequality v > 2/d, where d is the spatial dimension- 
ality of the system. In the opposite case, v < 2/d, the 
clean critical behavior is unstable, and the disorder ei- 
ther leads to a new, different universality class, or to an 
unconventional critical point, or even to the destruction 
of the phase transition. 

Another, less well understood consequence of quenched 
disorder is the formation of rare locally ordered regions 
in the disordered phase. For a transition occuring at a 
finite temperature, this can be explained in the follow- 
ing way. In general, disorder leads to the suppression of 
the critical temperature from its clean value T® to T c . 
In the temperature region between Tj? and T c the sys- 
tem does not show long-range order. However, there will 
be arbitrarily large regions which are devoid of impuri- 
ties and thus order locally. The probability of finding 
such regions usually decreases exponentially with their 
size, they represent non-perturbative degrees of freedom. 
These locally ordered regions are known as rare regions, 
and the order parameter fluctuations induced by them 
as local moments or instantons. Griffithso showed that 
the rare regions lead to a non-analytic free energy every- 
where in the temperature region between T® and T c , now 
called the Griffiths region or Griffiths phase. In generic 
classical systems this is a very weak effect, and the non- 
analyticity in the free energy is only an essential one. 
However, the Griffiths singularities become stronger if the 
disorder is spatially correlated. McCoy and WuB studied 
a two-dimensional Ising model where the disorder is per- 
fectly correlated in one spatial direction and uncorrelated 
in the other. In this model the rare regions lead to the 



divergence of the susceptibility at some temperature T x 
within the Griffiths region. 

A very interesting question is what is the influence of 
the rare regions on the critical behavior of a system. Dot- 
senko et al.Q studied this question for a weakly disordered 
classical ferromagnet. They found that the conventional 
theory of critical bchaviora in this system is unstable with 
respect to replica symmetry breaking. They also showed 
that the rare regions actually induce replica symmetry 
breaking perturbations and thus destabilize the conven- 
tional critical fixed point. While so far no final conclusion 
about the fate of the transition in the weakly disordered 
ferromagnet could be reached, the occurrence of replica 
symmetry breaking raises the possibility of an unconven- 
tional transition with activated scaling, as is believed to 
occur in the random field Ising madeB. 

For quantum phase transitional which occur at zero 
temperature as a function of some non-thermal control 
parameter, one expects an even stronger influence of the 
rare regions than for classical transitions. The reason is 
that a quantum model with uncorrelated quenched dis- 
order is effectively equivalent to a classical model with 
the disorder being perfectly correlated in one dimension 
(the imaginary time dimension) . Fishera investigated the 
critical behavior of a one-dimensional quantum Ising spin 
chain in a transverse field which is equivalent to the clas- 
sical McCoy- Wu model. He found that due to the rare 
regions the critical behavior is of the activated,form. This 
has beep-confirmed by numerical simulationsu which also 
suggestEj that this sort of behavior may not be restricted 
to one-dimensional systems. 

In two recent papersEJ we have considered the effect 
of rare regions on quantum phase transitions of itiner- 
ant electrons in d > 1. We have developed a systematic 
approach, representing the local moments by inhomoge- 
neous saddle point solutions of the field theory. The in- 
teraction between the local moments and the fluctuations 
leads to a new term in the effective action which is of the 
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form of annealed static disorder. In the case of the quan- 
tum antiferromagnetic transition this new term results in 
the destruction of the conventional critical fixed point if 
the number p of order parameter components is smaller 
than 4. No new fixed point could be identified, the sys- 
tem displays runaway flow to large disorder strength. On 
the other hand, for the quantum ferromagnetic transition 
the rare regions do not affect the critical behavior since 
a self-induced long-range interaction suppresses all fluc- 
tuations including those produced by the local moments. 

In this paper we apply the approach developed in Ref. 
|Tf| to a model of an itinerant antiferromagnet with an 
additional interaction term with cubic symmetry. This 
model is equivalent to a weakly disordered classical fer- 
romagnet with cubic anisotropy in which the disorder is 
perfectly correlated in some of the spatial dimensions but 
uncorrelated in the remaining dimensions. The conven- 
tional theory for this model (without taking rare regions 
into account) has ,bt;cn developed by Yamazaki, Holz, 
Ochiai and Fukudalij. 

The purpose for this work is threefold. We want inves- 
tigate (i) whether the conventional critical fixed point is 
stable under the influence of the rare regions. If it is un- 
stable we want to find out (ii) whether a new stable fixed 
fixed point exists which describes a rare region driven 
transition. Finally we want to study (iii) the influence 
of the rare regions on the fluctuation-driven first-order 
transition occurring in our system. The layout of the pa- 
per is as follows. In Sec. we derive the effective field 
theory by taking into account the disorder induced rare 
regions. In Sec. III. we carry out the renormalization 
group analysis. Finally, Sec. IV is left for a summary of 
our results. 



s G [4>] = \J dxdy 52M*)r(*-v)Mv) . ( 2 - lb ) 



£int [<t>] = U / dx ^ <j>j (x) <f)j (x) (x) (x) , (2.1c) 



5 cubi c[0] = A / dx X> 4 0*0 • ( 2 - ld ) 



Here we use a 4- vector notation to combine the real space 
coordinate x and imaginary time r, x = (x, r), J dx = 

J dx Jq^ T dr. The bare two point function, 

r(x-y,r-r') =r (x-y,T-T / ) 

+5(x-y)5(T-T')5t(x) , (2.2) 

consists of the deterministic part derived by HertzEl 
whose Fourier transform reads 



r (q,w„) = t + q 2 + \u n \ 



(2.3) 



and a disorder part in the form of a " random mass" term. 
Here q is the wave vector, oj n is a bosonic Matsubara 
frequency and 5t(x) is a random function of position and 
is endowed with the following statistical properties: 



(&(x))=0 , 



(<ft(x)<ft(y))=A£(x-y) 



(2.4a) 



(2.4b) 



II. AN EFFECTIVE ACTION FOR DISORDERED 
ANTIFERROMAGNETS WITH CUBIC 
ANISOTROPY 

A. The model 

In 1976 HertzB derived an order parameter field the- 
ory for the description of the antiferromagnetic quantum 
phase transition of itinerant electrons. Later this model 
was generalized to the dirty case by making the distance 
from the critical point a random function of positionEJ'E-X 
Here we consider an extension of this order parameter 
field theory by incorporating an additional (j) A term which 
possesses a (hyper-)cubic symmetry. 

In terms of the p-component order parameter field cf> 
(with components <pi) the total action can be written as 

S[<j>] =S G [4>] + S iut [4>]+S cuhic [<t>] , (2.1a) 

with the Gaussian part, the interaction part and the cu- 
bic anisotropic part given by 



B. Inhomogeneous saddle points and annealed 
disorder 

In the conventional approach, to critical behavior in 
systems with quenched disordera the disorder average is 
carried out at the [beginning of the calculation by means 
of the replica trickej. A subsequent perturbative analysis 
of the resulting, spatially homogeneous effective theory 
misses the rare regions we are interested in since they are 
non-perturbative degrees of freedom. 

We therefore follow the approach developed in Ref. [ll], 
and work with a particular realization of the disorder 
rather than integrating it out. Let us consider spatially 
inhomogeneous, but t ime -independent saddle point so- 
lutions of the action (|]l]) (time-dependent saddle-point 
solutions - if any - will always have a higher free en- 
ergy since the disorder is static). Depending on the sign 
of the cubic interaction term the structure of the saddle 
points in the p-dimensional order parameter space will 
be different. When A > the free energy is minimized 
by saddle point solutions that lie on the diagonals of a 
p-dimensional hypercube, while when A < the free en- 
ergy is minimized by solutions that lie on the axis of the 
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hypercubc. In cither case the modulus (f> sp of these min- 
imizing saddle point solutions fulfills the equation 

(t + 5t{x) - dl) |0 sp (x)| + 4u cS |0 sp (x)| 3 = , 

(2.5a) 



"off = 



u + I for A > 
u + A for A < 



(2.5b) 



Although sp (x) = is always a solution, there will be 
spatially inhomogeneousLSolutions if <5t(x) has sufficiently 
deep and wide troughaij. Let us now consider the Grif- 
fiths region, i.e. the region where the average distance to 
from the critical point is positive but where there are iso- 
lated islands which support a non-zero </> sp . If we have N 
such islands which are sufficiently apart from each other 
the global saddle point solutions may be written as 



N 

*£'>(x) = ${«}(x) = y>,(x) ffI 



i=i 



(2.6) 



where ipi(x) is a solution of ( |2.5| ) on the island / and crj 
is a unit vector in spin space (on one of the axis for A < 
or on one of the diagonals for A > 0). 

Since the direction of the order parameter on each of 
the TV islands can be chosen independently, (2.6) de- 
scribes an exponentially large number of degenerate sad- 
dle points, (2p) N for A < and {2P) N for A > 0. To be 
precise, the saddle points are not exactly degenerate due 
to the residual interaction of the (exponentially small) 
tails of the order parameter between the islands. The 
complicated structure of the free energy landscape con- 
nected with the existence of an exponentially large num- 
ber of almost degenerate saddle points will finally turn 
out to be responsible for the failure of the conventional 
approach. 

W e now consider fluctuations around the saddle points 
(2.6). Since the saddle points are separated by large free 
energy barriers an expansion around one of them will 
not give a good representation of the partition function 
of the entire system. Instead we will restrict ourselves 
to small fluctuations and simply add the contributions 
coming from all of the saddle points. Thus the partition 
function for a particular realization 5t(x) of the disorder 
can be written as 



Z[St(y 



E 



(2.7) 



Here f < indicates that the integration is restricted to 
small fluctuations ip only. 

We now carry out the sum over the saddle point con- 
figurations. The residual interaction between the islands 
will lead to slight deviations of the saddle point function 



from the ideal one given in (2.6). This is taken into ac- 
count by replacing the sum over the saddle points by an 
integral over a probability distribution 



P[<h\ 



a -i J dxC sp (<t>) 



(2.8) 



The temperature factor in the exponent reflects the fact 
that the_saddle points are classical (static) degrees of 
freedomtll. Expanding in powers of the fluctuations, we 
obtain the following effective action for the fluctuations 
ip (still for a particular disorder realization) 

S eS - S SP = S G [ V } + S int [(p} + S cuhic [cp] 

+ Tw I dxdy C(x,y)J2tf(x)rf(y) 



(2.9) 



higher order terms 



The correlation function C{x,y) measures, up to a con- 
stant factor determined by the precise form of £, whether 
x and y belong to the same island, and w — [(2 + 4/p)w,+ 
6A/p] is a positive constant. The w term is produced 
by the interaction of the fluctuations with the rare re- 
gions. It is our approximation of the effect of these non- 
perturbative degrees of freedom. Terms of higher than 
fourth order in ip also arise, but they are renormalization 
group irrelevant at both the Gaussian and the nontrivial 
fixed points of the conventional theory (see below) . 

Having identified the effects of the rare regions we now 
use the replica trickc3 to perform the quenched disorder 
average over St(x) which implies an average over position 
and size of the rare regions. The resulting effective action 
reads 

S eS [<p a (x)} = 

= ^EE f dxdyT (x-y)ip?(xM(y) 
+ U EE / dxdr(^(x,r)) 2 (^(x,r)) 2 
+ A EE / ^dr(^(x,r)) 4 
- A EE / d * dTdT ' (V?(x,r)) 2 (^(x,r')) 2 

a, (3 i,j 

- T *EE/ d * dTdT ' (V?(x,r)) 2 (^(x,r')) 2 



a,0 i,j 



(2.10) 



Here the first four terms are identical to the result of the 
conventional treatment. The 5th term has the form of 
static, annealed disorder and represents the interaction 
of the fluctuations with the rare regions in the Griffiths 
phase. For more details of this derivation see Ref. [Tl]. 
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III. RENORMALIZATION GROUP ANALYSIS 



B. Fixed points and their stability 



A. Flow equations 



We first consider the effective action (2.1C) at tree 
level. As usual, let us define the scale dimension of a 
length L to be [L] = — 1, and that of imaginary time r 
to be [r] = — z with z being the dynamical critical expo- 
nent. We first analyze the Gauss ian fi xed point. From 
the Gaussian part of the action ( 2.1 0| ) we see that ui n 
scales as q 2 , implying that 2 = 2. The scale dimension 
of the field is [<p] = d/2. Power counting for the inter- 
action and disorder terms of the action gives the scale 
dimensions of it, A, A and w as [u] = [A] = [w] = 2 — d 
and [A] = 4 — d. Here we have used the fact that in 
Matsubara formalism the temperature scales like a fre- 
quency, [T] — z. Consequently, u, A and w are irrelevant 
for d > 2, while A is irrelevant only for d > 4. This im- 
plies that in the physical dimension d = 3 the Gaussian 
fixed point is unstable, and we must carry out a loop 
expansion of the effective action (2.10) close to d — 4. 
All terms of hi gher o rder in (p that arise in addition to 
those given in (2.10) have negative scale dimensions at 
and close to d — 4. Thus, they are irrelevant by power 
counting with respect to both the Gaussian and the con- 
ventional non-trivial fixed points^,,— ■ ■— ■ 

As in the conventional theoryli3o't3 we carry out the 
perturbation theory in d = 4 — e spatial dimensions and 
e T time dimensions. In this way the perturbation ex- 
pansion becomes a double expansion in terms of e and 
e T . The renormalization group flow equations are ob- 
tained by performing a frequency momentum shell RG 
procedure.Eij To one-loop order, we obtain the following 
flow equations, 



du 

~dl 

dX 

~dl 
dA 

~dJ 

dw 

~df 



eu - 4(p + 8)u 2 + 48uA - 24mA, 
iX - 36A 2 + 48AA - 48wA, 



(3.1a) 
(3.1b) 



eA + 32A 2 - 8(p + 2)uA + 8pAw - 24AA, (3.1c) 
ew + Apw 2 - 8(p + 2)uw + 48Aw - 24Aw. (3. Id) 



Here we have defined e = e — 2e T . Of course, also the 
distance t from the critical point will be renormalized. 
However, we only consider the flow on the critical sur- 
face t = since we are interested in the stability of the 
critical fixed points. Note that the coefficient of the rare 
region term w only couples to A. The flow of u and A is 
only indirectly influenced by the rare regions (via a mod- 
ification of the flow of A). This will be important later 
on. 



The flow equations (3.1) possess sixteen fixed points. 
Their fixed point values are given in Table |[ the eigenval- 
ues of the corresponding linearized renormalization group 
transformations are listed in Table || For eight of the 
sixteen fixed points (Nos. 1-8 in Table ||) the fixed point 
value of the rare region term is w* = 0. These fixed 
points have already been studied in Ref. ^ using the 
conventional approach. In the following, we concentrate 
on the case e > and e = e — 2e T < relevant for the 
itinerant quantum antiferromagnet. 

We first consider the dirty Heisenberg fixed point (No. 
6) and the dirty cubic fixed point (No. 8). These are the 
stable fixed points of the conventional theory for the cases 
of p < 4 and p > 4, respectively. Analyzing the stability 
matrix for the dirty Heisenberg fixed point shows that it 
is unstable since the eigenvalue is positive for p < 4. 
In contrast, the dirty cubic fixed point remains stable 
for p > 4 since all eigenvalues of the stability matrix are 
negative. Thus we conclude that the rare regions destroy 
the conventional dirty Heisenberg critical behavior for 
p < 4 while they do not influence the conventional dirty 
cubic critical behavior for p > 4. 

We now turn to the new fixed points with w* ^ (Nos. 
9 - 16 in Table |). Fixed points 9, 11, 13 and 15 are un- 
physical because their fixed point values w* are neg ative . 
Since the bare w is positive and according to eq. ( 3. Id ) 
the flow cannot cross the (w = 0)-plane these fixed points 
can never be reached. Depending on the number p of 
order parameter components the remaining fixed points 
(Nos. 10, 12, 14, and 16) are either also unphysical, or 
they are unstable. Consequently, for p < 4 and to one- 
loop order there is no stable fixed point. Renormalization 
group trajectories which in the conventional theory would 
go to the dirty Heisenberg fixed point show runaway flow 
to large disorder strength. This runaway flow could ei- 
ther indicate a unconventional phase transition, e.g. an 
infinite disorder critical point as in the one-dimensional 
random Ising modelEl or a percolative rather than a ho- 
mogeneous transition or even a destruction of the phase 
transition. Within the present approach we cannot be 
decide between these alternatives. 

The influence of the rare regions on the stability of 
the fixed points in our model is similar to that in the 
isotropic caseO. For p < 4 the conventional fixed point 
is destroyed in both models. For p > 4 the conventional 
fixed point is stable. In our model this is the dirty cubic 
fixed point while in the isotropic case this stable fixed 
point is the dirty Heisenberg fixed point. 



C. The fluctuation-driven first-order transition 



In addition to the continuous phase transitions associ- 
ated with the critical points discussed above there is also 
the possibility for a first-order transition in the model 
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No. 










FP values 








11* 




A* 


A* 


w* 


1 


















2 




e/4(p + 8) 













3 









e/36 








4 




e/12p 




?(p - 4) /36p 








5 












-e/32 





6 


(3e 


- 2e)/16(p - 


1) 





[(p + 8)e-2(p + 2)e|/64(p-l) 





7 









C(6 1/2 ) 


e>(6 1/2 ) 





8 


(36 


- 2e)/24(p - 


2) 


[(3e-2e)(p-4)]/72(p-2) 


[3pe - 4(p - l)e]/96(p - 2) 





9 















-e/4p 


10 




e/4(p + 8) 










[(p-4)e]/4p(p + 8) 


11 









6/36 





-e/12p 


12 




c/12p 




[(p - 4)e]/36p 





(p-4)67l2p 2 






u 




U 




^ze — o€) / lop 


14 


(36 


- 2e)/8(10 - 


P) 





[O + 8)e-12e]/32(10-p) 


[(3e-2e)(p-4)]/8p(10-p) 


15 









(3e - 2e)/72 


(9e- 12e)/288 


-3(36- 2?)/72p 


16 




(3e - 2e)/48 




(3e-2e)(p-4)/144 


(3pe - 2(p + 2)e)/192 


(3e-2e)(p-4)/48p 



TABLE I. Fixed points of the flow equations, p is the number of order parameter components. 



considered here. Let us first discuss the mechanism for a 
clean system and discuss the effects of disorder and rare 
regions later. 

According to a m ean-field stability analysis of the ef- 
fective action ( 2.10| ) with A = w = the inequalities 
u + A > (for u > 0) and u + X/p > (for u < 0) have 
to be fulfilled for the theory to be stable. Now consider 
a bare theory with « < 0, A > or u > 0, A<0 but still 
fulfilling the above s tab ility conditions. In these cases 
the flow equations ( |3.l| ) can lead the renormalization 
group trajectories to the mean-field unstable region. 1L, 
indicates a fluctuatina,-driven first-order transitiorM 
It was later shownEJ'EJ that the fluctuation-driven first- 
order in this model survives the presence of quenched 
disorder, at least within the conventional theory. Let us 
now consider the influence of the rare regions. As al- 
ready mentioned, the rare region coefficient w does not 
couple into the flow equations for u and A but only into 
the flow equation for A. Thus a renormalization group 
trajectory going to the mean-field unstable region within 
the conventional theory will generically also do so in the 
presence of rare regions, the only modification being a 
different disorder value at the stability boundary. 

Therefore, we conclude that the fluctuation-driven 
first-order transition also occurs when taking the rare re- 
gions into account. However, since the rare regions mod- 
ify the flow of the disorder strength A, the boundaries of 
the first-order region may change. 



IV. SUMMARY AND CONCLUSIONS 

We have investigated the influence of rare regions on 
the quantum phase transition of a disordered itinerant 
antiferromagnet with cubic anisotropy. The local mag- 
netic moments forming on the rare regions in the Grif- 
fiths phase generate a new term in the order parameter 
field theory which has the form of static annealed disor- 
der. We have found that for order parameter dimension 
p > 4 this new term does not change the critical be- 
havior, which is characterized by the dirty cubic fixed 
point. In contrast, for p < 4 the rare region term renders 
the conventional critical fixed point unstable. The renor- 
malization group trajectories show runaway flow to large 
disorder. Within our approach which is essentially per- 
turbative, even though it includes some non-perturbative 
degrees of freedom (the local moments) we cannot deter- 
mine the ultimate fate of the transition. It could be an 
unconventional phase transition, e.g. an infinite disorder 
critical point or a percolative rather than a homogeneous 
transition or even the destruction of the phase transition. 
We have also found that the fluctuation-driven first-order 
transition occurring in this model remains qualitatively 
unchanged by the rare regions, while the precise position 
of the first-order region in parameter space will change. 

The authors acknowledge helpful discussions with D. 
Belitz, J. Cardy, and T.R. Kirkpatrick. R.N. thanks the 
hospitality of T U Chemnitz during two vistits where part 
of the research was performed. This work was supported 
in part by the DFG under grant nos. SFB393/C2 and 
Vo659/2, by the NSF under grant no. DMR-98-70597, 
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eigenvalues 

ei e 2 es e 4 



e e e e 

-g (p-4)e/(p + 8) e-2(p + 2)g/(p+8) -(p-4)g/(p + * 

e/3 -g e-2g/3 e/3 

-g -g(p - 4)/3p e - 4g(p - 1) /3p -e(p - 4)/3p 



eigenvalues not calculated since FP is unphysical 

-a + £^ ~ A -ft* (p-4)(3e-2g)/4(p-l) -(p - 4)(3e - 2g)/4(p - 1) 

~ E tS^ -(3 £ - 2g)(p - 4)/6(p - 2) -(3. - 2e)(p - 4)/6(p - 2) 



eigenvalues not calculated since FP is unphysical 

-g (p_ 4 )g/(p + 8) e-12g/(p + 8) (p-4)g/(p + 8) 

eigenvalues not calculated since FP is unphysical 

-e -e(p-4)/3p e-2g(p + 2)/3p g(p - 4)/3p 



eigenvalues not calculated since FP is unphysical 

"^uST 15 " C 4"(uST° (p-4)(3e-2g)/2(10-p) (p - 4)(3e - 2g)/2(10 - p) 

eigenvalues not calculated since FP is unphysical 

-G+^G 2 -H -G-y/G^-H 



_>-l 



(3e-2g)(p-4)/12 -(3e-2g)(p-4)/12 



TABLE II. Eigenvalues of the corresponding linearized RG transformation, p is the number of order parameter components. 
A, B, C, and D are defined as A = (p + 8)e - 2(p - 4)g, S = 16(p - 1) (3e - 2g) [(p + 8)e - 2(p + 2)g], C = (p + 8)e - 2(p - 4)g, 
D = 8(10-p) (3e-2g) [8e-12g+pe]. Analogously, E = 3pe + 2(p-4)g, F = 24(p-2) (3e-2g) [4g+3pe-4pg], G = 8g+3pe-2pg, 
H = 48(3e - 2g) [-4g + 3pe - 2pg]. 



and by EPSRC under grant no. GR/M 04426. 
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